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SUMMARY 


The present study evaluates a higher-order modal method proposed by Leung 
for transient structural analysis entitled the force- derivative method. This 
method repeatedly integrates by parts with respect to time the convolution- 
integral form of the structural response to produce successively better approxi- 
mations to the contribution of the higher modes which are neglected in the modal 
summation. Comparisons are made of the force- derivative, the mode- displacement, 
and the mode- acceleration methods for several numerical example problems for 
various times, levels of damping, and forcing functions. The example problems 
include a tip-loaded cantilevered beam and a simply-supported multispan beam. 

The force-derivative method is shown to converge to an accurate solution in 
fewer modes than either the mode- displacement or the mode-acceleration methods. 
In addition, for problems In which there are a large number of closely- spaced 
frequencies whose mode shapes have a negligible contribution to the response, 
the force derivative method is very effective in representing the effect of the 
Important, but otherwise neglected, higher modes. 
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cross-sectional area 
damping matrix 
modulus of elasticity 
error norm (see eq. (17)) 
i th modal forcing function 
moment of inertia 
stiffness matrix 
beam length 
mass matrix 
bending moment 

number of modes used in truncated modal summation 
total number of degrees of freedom 

applied non-dimensional force vector (Q=qL 2 /EI or Q=qL 3 /EI) 
applied load 
distributed load 
frequency ratio 

u>i 

shear force 
time 

normalized time T = u> 0 t 

structural displacement response vector 

cartesian coordinates 

i th modal displacement function 


Greek : 

Y 

{«} 

Ci 

y(t) 

P 

{♦} i 
[*] 
[n] 

“df 

u>i 


Subscripts : 

c 

m 

o 


forcing frequency 

vector of difference between converged response and approximated 

response using m modes 

i th modal viscous damping factor 

unit step function 

mass density 

dummy variable of integration 
i th natural vibration mode 
modal matrix 


diag [ux ,u> 2 • • diagonal matrix of circular natural frequencies 
itt> circular frequency of the damped free vibration 
i th circular natural frequency 

pAL^ 


normalizing frequency u> 0 * 



converged 

response using a truncated set of modes (1 to m) 
i ni ti al 


Superscripts: 

(•) 

(n) 


differentiation with respect to time 
n th derivative with respect to time 



INTRODUCTION 
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Dynamic analysis of complicated structures which are modeled as discrete 
multi degree-of- freedom systems often requires the solution of very large systems 
of equations. Reducing the order of such systems is highly desirable from the 
standpoint of increased computational efficiency. Some of the many methods for 
reducing the order of discrete multi degree- of- freedom dynamic systems include 
mass condensation methods (refs. 1-2) and reduced basis methods (refs. 3-7). 

The reduced basis methods use the expansion theorem and a truncated set of basis 
vectors (e.g. undamped free- vibration modes (modal methods) or Ritz-vectors) to 
approximate the dynamic response. Some problems such as large space structures 
require a large number of basis vectors to accurately represent the dynamic 
response and in some cases, such as when singularities occur in the loading, 
convergence of a solution is not guaranteed. Also, for problems which require 
the derivatives of the response with respect to a design parameter (e.g. system 
identification or optimization problems) not only is the size of the problem 
increased but the convergence of the derivative equation is not generally 
guaranteed. 

When a reduced basis method uses the natural vibration modes of the 
structure the method is referred to as a modal method. Two of the most 
widely-used modal methods are the mode- displacement method (MOM) and the 
mode- acceleration method (MAM). Comparisons of the MDM and MAM methods (ref. 6) 
indicate that the MAM converges to an accurate solution with fewer modes than 
the MDM. One reason for this improved convergence is that the MAM incorporates 
a pseudo- static response (the inverse of the stiffness matrix multiplied by the 
applied forcing function) which approximates, to some degree, the flexibility of 
the higher modes which are neglected in the modal summation. The work of 
reference 7 reveals that the MAM method can be derived by integrate ng-by-parts 
with respect to time the convolution- integral form of the original MDM method. 

It is also shown in reference 7 that higher-order modal methods may be obtained 
by further integration- by- parts. Thus, integrating the convolution integral two 
more times produces a higher-order modal method (ref. 7) than the MAM which is 
called here the force-derivative method (FDM). It is called the 
force- derivative method because, analogous to the MAM, the FDM produces a term 
which is a function of the forcing function and additional terms which are 
functions of the time-derivatives of the forcing function. These additional 
terms produce successively higher-order approximations to the higher, neglected 
modes. It is shown that for the case of zero damping the MAM is a higher-order 
method than the MDM, and can be derived by integrating the convolution integral 
by parts two times. The purpose of the present study is to evaluate the use of 
higher-order modal methods which have the ability to significantly reduce the 
number of degrees- of- freedom (modes) necessary to accurately represent the 
transient structural response. The present study extends the work of reference 
7 to include modal damping and clarifies the relation of the MDM and MAM to the 
FDM. 


The rate of convergence or accuracy of each of the modal summation methods 
(MDM, MAM, and FDM) is investigated with respect to the differentiability of the 
forcing function, the level of damping, and the time at which the response is 
calculated. A relative, spatial error norm is used to measure accuracy and con- 
vergence of the transient response. The forcing functions were chosen to illus- 
trate the effect of a continuous forcing function (with respect to time) which 
has continuous higher derivatives and a discontinuous forcing function (a unit 
step function) on the solution. A quintic function of time was selected to 
illustrate what happens when the force or one of its derivatives vanish at some 
point in time. The convergence for a sinusoidal forcing function is also 
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studied. A series of numerical examples has been selected to illustrate the 
adequacy and/or inadequacy of each method: (1) a uniform cantilevered beam sub- 

ject to a tip loading condition of (a) a unit step load input, and (b) a quintic 
function of time; and (2) a simply- supported multi span beam (10 equal length 
spans) subjected to (a) a uniformly distributed load varying as a quintic 
function of time and b) two concentrated loads located about the center of the 
first span and varying as a quintic function of time. 

UNIFIED DERIVATION OF MODAL METHODS 

The equations of motion, in matrix form, of an n-degree of freedom system, 
together with the initial conditions, are given by 


[M]{U( t)} + [C]{U(t)} + [K]{U( t)} = { Q ( t)} (1) 

{U(o)) = fu 0 ), {U(o)| . (U 0 ) 

where [M], [C], and [K] are the mass, damping, and stiffness matrices of the 
structure; U(t) and Q(t) are the displacement and force vectors, respectively, 
and a dot represents differentiation with respect to time. For certain forms of 
[C], equation (1) can be transformed into a set of uncoupled equations of motion 
using the undamped normal modes, {^K, of the structure and the expansion 
theorem, namely, by 

{u( t)} “jX UMi(t) (2) 


where the natural vibration modes are obtained from the solution of the 
following eigenvalue problem 

CK3{4»} i = 0)^[M]{*} i (3) 

The modes may be normalized such that 

{♦}|[M]{(|.} i = 1. and then { 4 } J[K]{<|>} w* (4a) 

Using the modal matrix [*] which contains all n modes equations 4(a) can be 
written as 

[*] T [M][*] = [I] and 0] T [K][*] = [n 2 ] (4b) 

In modal coordinates, equation (1) can be expressed 
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V*> + ♦ »*Y f (t) ■ F,(t) 1=l,2,3...n < 5 » 

Y 0j * t*l>ltu 0 }. » 0j * !♦)>]{ ij 

where the i th modal force Fj(t), Is defined by 

F.(t) = {*}] { Q ( t) } ( 6 ) 

and 5 , is the i^* modal viscous damping factor. 

The mode-displacement method (MDM) calculates m- values of Y-j(t) by solv- 
ing the first m-equations (eq. 5) and substituting this into a truncated form of 
equation ( 2 ). 

m 

{ U( t)} = J {♦} - ( t) where m£n ( 7 ) 


The solution of equation (5) can be written by using a convolution integral 
as f ol 1 ows 


-CjOjjt i 

Y,(t) = e [cos U dj t> Y 0j+ — (Y 0) + c,»,Y 0 j )s1n 

I t “5 ( t-T ) 

+ -jjj— / e sin (tD^ ( t-x )){F^(x)}dx 


(8) 


d. o 


where 


w d. = “i^ 1 ‘ 


Substituting equation ( 8 ) into equation (7) gives 
tu(t)| -.!.{{♦!, te <cos u d t) Y + 1 (Y + c,»,Y 0 .)*in (.„«)] 

1 1 a. 1 1 l 


(9) 


i=l 


i 

1 t -?iO>j(t-T) 

+ {*} i [ ^/e 11 sin (a, d (t-x)) F^xJdx]} 

w dj o u i 


Equation (9) is an equivalent representation of the MDM. Integrating the convo- 
lution integral of equation (9) by parts twice with respect to time and assuming 
the forcing function and its first derivative are continuous over the interval 
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m 


-w 


i 


|u(t)) -H {♦(,[• ' ■ (cos u d t) Y 0 . +;My ♦ CiVo.'s'" 
1 = 1 


1 1 w i i 1 


{♦} 


i [e’ Ci “ it ^7 (^- sin (a» d t) + cos U d t))]{<|>}J{Q 0 } (10) 


w a 

i d. 


_ u t 2 J. 

+ {♦} t Ce Cl 1 1 - d j ^ sin U d t) + Zc.-uj.cos (a> d t) )]{♦} J{Q 0 } 

f a) j i 1 

w t d i 

? 2 

. t -C.U.(t-T) /C.-W4Y _ “d. 

+ {♦} -ft — 2T / e H iT ^ sin 


2*; ..a^COS ((Oj ( *t — t ))}{♦} ^{Q^(t )} dr ] 


+ - {4>} ^ C — } 

Ulj ID 4 


2c 1“i 


i Tr A, 


If the forcing function or its first derivative were discontinuous at some point 
in time the appropriate jump conditions would have to be Included in 
equation (10). 

If all the modes are used in calculating the last two terms of equation 
(10) these terms can be represented as 


f {♦} i (^){+}T{ Q (t ) } = [*][tf 2 ][*] T lQ(t)} =[*][[*] T [K][*]] _1 [*] T = CK3 _1 {Q(t)} 

(11) 


1=1 ' o >i 


and similarly 


I {♦} 1 (^X-M+ i } T {Q(t)} - C*][^-][*] T {0(t)} = CK3" 1 Cc3CK] _1 {Q(t)} 

1=1 0)^ U)^ 


where [*] is the complete (nxn) matrix of eigenvectors and [q“ 2] is the 


7 


2c .u. 


11 1 ^ ii 

diagonal matrix -=y » • and [ — jp-] is the diagonal matrix 

U )^ 1&2 < a ) 4 


2^ 2 52 u, 2 _ 2c„» n 

[ — r — j- — r !• 

<l)i 0>o <»>» 


The term, [K]‘l{Q(t)}, is called the pseudo-static response in reference 6. 
Hence, if all the modes are used in the last term of equation (10) it is 
expressed as 

m “5 i "t i • 

{u(t)} (cos (w d t) Y + — (Y 0 + c i u. i Y 0 )sin U d t))] 

1 1 a. 1 1 1 


1=1 


‘1 


+ 7? sin ( U(J t) ♦ cos (« d .t))](*}|(Q 0 ) (12) 

i dj 1 i 

+ {♦jjte ** 1 1 ~ d i iSin (w d t) + Zc^^cos (w d t))]{+}{{Q 0 } 

T" Wj a i U 1 


2 w 2 

, t ”C jWj( t-x ) iZi<a*\ 

+ e 1 1 {(. U 0 -- - -- l )si n (w d (t-x 

<a»j 0 1 


)) 


+ 2c.gw.-cos (« d (t-x))}{4»}^{Q^(T)}dx] } 


+ [k3 _1 {q( t)} - urkciw^Qtt)} 

For the undamped case, equation (12) can be easily shown to be equivalent to the 
HAM given as 


{U(t)} * [K] _1 {Q(t)} -l {^} i [(-i 2 -)Y*(t)] (13) 

1=1 U) g 

As shown by equation (13) the MAM uses the modal expansion to approximate the 
deviation from a pseudo steady-state response. It can be expected to perform 
better than the MDM for low frequency excitation, where the response can track 
the load well. The pseudo- static term, by making use of the Inverse of the 
stiffness matrix of the structure, is a first-order approximation to the flexi- 
bility of the higher, neglected modes. 
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The force derivative method makes use of the fact that successive Integra- 
tions by parts of the convolution Integral produces terms which are functions of 
the natural frequency to successively higher negative powers and terms which. If 
all the modes are used In the summation, can be represented by the system mass, 
stiffness, and damping matrices and the forcing function and Its derivatives. 
Hence, higher-order methods can be developed by successively integrating the 
convolution Integral by parts. 


Integrating the convolution Integral of equation (10) two more times, one 
obtains 

m i • 

{U(t)} = {♦} i [e {cos u d t) + — (Y 0 + c^Y^sin U d t)}] 


1=1 


*1 


d 1 1 


'1 


'1 


+ Ulfte sin U d t) + cos (u> d t) ){*}}{ Q q }] 

dj i i 


(14) 


“CjiDit « (<i)jC{) - d. _ . 

+ {+}^[e — ^(( ) si n U d t) + 2c.fCD.jCos (cD d t))]{4>} ^{Q 0 } 

W 1 w d i i i 

2 

" C i 0 , 1 t 1 O T • • 

+ [*)£* ~ > s1n ( “d + ^- 4 4 )cos S 

u.| d^ 1 1 

-C^u.t . l-8c? + 8;} (8c?-4cJ 

+ HhCe 1 {-t-(( — ) si n («, t) + — 2 — 

11 1 ‘ 4 ia A d., <d 4 

“1 


“ ~~ cos U di t))}]{ + }{{Q 0 } (3) 


+ Mi a 


t e -?i“i (t - T) (l-8c}+8c}) 


CD, 


“d, 


sin (o»j ( t — t ) ) 
a 1 


(8^-4;,) T m 

+ — * - cos U d (t-T))]{+}J{Q(T)} (4, dr] 


+ {♦} i C-^3Hl{{Q(t)} - {♦} 1 [-4]{*}{{Q(t)} 

CD .j CD ^ 

(8c?-4c-) T # 3 \ 

- ■] {♦}{ {Q(t)} 3 } 


a-4c?) t •• 
{♦ lit — r -]{ 4 } I { Q ( t )} 
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where 

(3) (3) 

{Q(oM - (Q 0 ). f Q Co)} ■ (Q 0 |. )Q(o)| • }Q 0 } » and (Q(o)} - }Q 0 } 

If all the inodes are used In the modal summation of the last four terms they can 
be represented as 

l s t K l 1 

1=1 U>. 


1-1 U) • 


>1 

n (1-45^ ) x i i 1 _l 1 

I Mil r-]{*}I * [K] 1 - [K] *[C][K] nCICK] 1 

1 = 1 Cl>^ 


-1 


(15) 


and 


(8^4) 

Uli [.iipi ] U1 T = CK3" 1 [C]CK}- 1 [C3CKl _1 CC3CK] _l - CK] _l tC3CK] _l tM3[K] _1 

“1 

Using equation (15), equat1on(14) becomes 

in —5 jt 1 • 

{U(t)} -j { {♦} i Ce {cos (u> d t) + — (Y Q +c 1 « i Y 0 )sin U d t)}] 
i=l 1 d d 1 1 1 


+ {♦} -f te Cia>i t - 1 { ; p sin (u, d t) + cos (« d t)}{*}{{Q 0 }3 
a>^ 1 i 


(16) 


{♦life — )s1n U d t) + 2c cos U> d t)} ]{♦}{{ Q q } 

“1 w d. 1 1 


1 

.2y 


“Cjto.t 1 c ,*<0 * (3-4£ . ) o T •• 

+ {♦} i Ce 1 1 {-^(-2-L — )s1n U d t) + (l-4 c 2 )cos («o d t))}]{ + }]{Q 0 } 

d,. 1 1 


IDj 
d 1 

- 5 iV, 1 ,, 1 - 8 4 + 


(8t 8 -4 Ci ) 


|T,„ ,(3) 


{♦)<!> ' ' (-!*<(— Main U d .t) t—i-J-cos U d t) )}]{♦) j(Q 0 ) 

ci) ^ d* 1 i 1 



10 


t / c l“l (t_T ) (l-8c?+8c^) 

+ {♦}*£/ $ C jr sin Ki. (t ' T 

0 u> ^ 


V 


)) 


(85.-4tj} T /*\ 1 

+ • ~ •— cos U„ (t-T»]{ + )I{Q(T)) 1,, <h] ) + CKl -1 {0Ct»l 

w i 1 


- [KrHcDtKr^Qtt)} - [ [K] _1 [M][K] _1 - CK3" 1 CC][K]' 1 CC3CK3" 1 3{Q(t)} 

- [ CK3” 1 CC3[K]" 1 CC3CK]“ 1 CC3[K]- 1 - [K]‘ 1 [C3[K]" 1 [M][K]“ 1 ]{Q(t)} (3) 


Hence, analogous to the MAM for zero damping which retains a first order approx- 
imation to the effect of the neglected modes by the [K]" 1 Q(t) or 
UHh ][*] T {Q(t)} term, the FDM offers a second-order approximation to the 

neglected modes by including the [K]" 1 [M][K]’ 1 {Q (t)} or U][n~*]i>] T {Q (t)} 
term. Results of reference 7 indicate that this higher-order- approximation can 
significantly reduce the number of modes required for convergence, depending on 
the nature of the forcing function. Higher-order approximations than those pre- 
sented for the FDM (using four integrations by parts) can be developed by suc- 
cesively integrating the convolution integral by parts several more times. 

IMPLEMENTATION 

The above equations were programmed for beam examples using analytical 
expressions for the mode shapes {<!>}.. and modal displacements and their deri- 
vatives. Modal vectors assume 51 equally- spaced points along the length of 
beam. Moment and shear forces were calculated from M(t)=El3 2 U(t)/3x 2 and 
S( t) a El3^u( t)/3x 3 . For most cases a converged solution was assumed to 
occur using 30 inodes and, hence, a total of 30 modes were used to approximate 
expressions such as 

do -1 - )!♦){, ocrkwcur 1 = 3 ° 

i=l ( 1 ).. i=l o>,j 


Whenever necessary a total of 50 modes is used in the above expressions to 
represent a converged solution. 

The following section investigates the effect of various forcing functions, 
load distributions, and damping levels on the transient response of a uniform 
cross-section cantilevered beam, and a uniform simply- supported multispan beam. 



Cantilevered Beam With Tip Loading 
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The first problem studied is a uniform cantilevered beam under tip loading. 
The first 30 natural frequencies of the beam are listed in table 1. The beam is 
subjected to various levels of modal damping (same for all modes) and a variety 
of loading conditions which are described below. 

Case 1: Q(T)=100Q(T 4 -T 5 ) .- This problem was presented in reference 7 for the 

case of zero damping Ui=0). The function or one of its derivatives vanishes 

• • • 

at yarious times: Q(T)=0 when T=1.0, Q(T)=0 when T=0. 8, Q (T)O when TO. 6 

q(3)(t)o when TO. 4, and q(4)(t)o when TO. 2 and this fact affects the 
convergence of the method as will be shown subsequently. As shown in figure 2 
for c i=0.05 and T=1.2, the FDM offers an improvement in accuracy of several 
orders of magnitude in the error norm over either the MDM or the MAM methods. 

The advantage in using higher-order modal methods (either the MAM or the FDM) 
lies in the ability of those methods to approximate the flexibility of the 
higher, but neglected, modes with terms which are functions of the stiffness, 
mass, and damping matrices and the forcing function and, in the case of FDM, its 
derivatives (see for example eqns. 14 and 16). The FDM offers a higher-order 
approximation to the neglected modes by the use of additional terms in addition 
to the pseudo-static response ([K]" l {Q(t)} ). Figure 3 is a plot of the moment 
error norm (using 5 modes) as a function of time; at T=1.0, the moment error 
norm associated with the MDM is equivalent to the MAM value because at time 
T=1.0, Q(T)=0 and there is no difference between MDM and MAM (see eq. 12). A 

similar dip in the MAM error occurs at time T=0.8 which corresponds to a time 
• 

when Q(T)=0. These narrow regions where there is a sharp increase in solution 
accuracy can be anticipated a priori from the zeroes of the forcing function and 
its derivatives. 

A comparison of displacement, moment, and shear errors for 4-j=0. and 
times T=0.6 and 1.0 are shown in figures 4a and 4b, respectively. The error 
associated with the displacements is lowest, the moment errors are greater and 
the shear errors are the largest. This occurs because the moments and shears 
are functions of successively higher spatial derivatives of the displacements 
(ref. 6). When T=0.6 and 4^=0., the MAM and FDM methods coincide (fig. 4a) 

• • 

because Q (T)=0. When T=1.0 the MDM and MAM methods coincide (fig. 4b). From 
figure 4b, when T=l. and 4i=0 the only difference between the MDM and MAM 

methods and the FDM lies in the term {$} 4 ] Q(T)} (eq. 16). 

“1 

The difference in errors between the methods is less at T=1.0 (Q(T)=0) than at 
T=0.6 (fig 4a) where the difference between the methods lies in the 

gH+J i{Q(T)} term. The higher-order terms are functions of the 

“1 

frequencies to successively higher negative exponents and, hence, should have a 
negligible effect as higher modes are used providing the time- function term does 
not grow proportionally. 



12 


RESULTS AND DISCUSSION 


The convergence of each method (number of modes versus accuracy of the 
transient response) is measured by using a relative, spatial error norm. The 
error norm, e, of an approximation to the displacement vector based on the first 
m modes {U} m is given as 


V{u} T J 


1«L 


(17) 


where 


{6} = {U} c - {U} m 

{ U} m and {U} c represents a converged solution (for most of this study 
the first 30 modes are assumed to provide a converged solution). Similarly, 
moment and shear error norms are calculated by substituting the moment or shear 
response vectors in equation 17. 

The effectiveness of the error norm, e, as a measure of accuracy and con- 
vergence in quantifying the global error associated with each of the modal 
methods is demonstrated for a cantilever beam problem. For this problem 
b ) q = /EI/pAL 4 and time t is normalized such that T = w 0 t. Figures la and 
lb demonstrate the use of the error norm in describing the accuracy of 
displacement and moment distributions for a tip load (Q(T) = 100Q(T 4 -T 5 )) at a 
normalized time T = 0.4 and for ci = 0.05. In figure la only one mode is 
used, and for the MDM e = 0.045 and displacement errors are noticeable, and MAM 
and FDM results having error norms of e = 0.00627 and e = 0.0000317, 
respectively, are indistinguishable from the converged solution (MDM( 30 ) ) . 
Similarly, for the normalized moment distribution (fig. lb) the MDM (using the 
first two modes) has an error e = .128, the MAM and FDM (both using only the 
first mode) have errors of e = .0379 and e = .00019 respectively. There is a 
qualitative improvement in the solution (response distribution) as the error 
norm e decreases in magnitude. 

The rate of convergence of each of the methods (MDM, MAM, and FDM) is 
expected to depend on the nature of the forcing function, the level of damping 
and the time at which the response is calculated. The forcing fuctions were 
selected to investigate the effects of continuous forcing functions with vanish- 
ing higher derivatives at various times ( Q(T) = 1000(T 4 -T 5 ) ), the effect of a 
discontinuous forcing function, a unit step at t = 0 (Q(T) = u(T)), and the 
effects of a sinusoidal forcing function (Q(T)=sin yT). It is assumed that for 
all forcing functions Q(T) = 0 for T<0. Several example problems, described 
below, were selected to evaluate the accuracy of each method. 


/ 
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The effect of damping on the accuracy of the response is shown in figures 
5a and 5b. Increasing the modal damping does not always Increase the 
accuracy of the MAM as suggested in reference 6. For the case studied in refer 
ence 6, a uniformly loaded cantilevered beam subjected to a step loading, the 
accuracy of the MAM is enhanced in the presence of damping as can be seen from 
equation (16). Since all the derivatives of the forcing function vanish, the 
only terms remaining are the pseudo-static response and a term which is a func- 
tion of 0- 5 l“it. Hence, as ^ increases, the relative importance 
of this term as compared to the pseudo-static term ([K] _1 {Q(t)} ) decreases and 
therefore the accuracy of the MAM increases. 


For other forcing functions, terms such as {<p} 


r 2 ** i“ i 

i L r 

0 ) • 


■] l»){|Q(t)| 


increase 


in importance as Cf increases and so the effect of ?i on accuracy is 
more complex. For example, the case of a linearly varying forcing function, 

Q ( t) = t, the damped system response can be represented, assuming zero initial 
conditions, as (see eq. (16)): 


m 


{U(t)} * l {♦} i [e 1 ' ±- 4 [( 

1 = 1 U) j 


' ? i“l* 1 r / (u> i c 1 )2 " U>d i 


•)sinU d t)+ 2c i u i cos(a) d t)] { 4 »}^{Q q } ] 


J i 


T * 

1 1 I (ni 


♦ [♦][a- 2 ][*l T (Q(t» - [•H-H-H*]' {Q(t» 

“1 


From this equation, the MAM assumes all the modes are used in the second to last 
term and it is represented as the pseudo-static response [K] -1 {Q(t)}. The 
remaining terms are represented by a truncated modal summation. The first term 
decreases exponentially as damping (tj) Increases but the last term increas- 
es proportionally to C-j. If the magnitude of this last term does not 
decrease with respect to the pseudo-static term the increase in damping will not 
necessarily result in an increase In accuracy as in the unit step function case. 

Case 2: Q(T) = u(T) .- This forcing function is discontinuous at time T=0 and 

hence the integration by parts of the convolution integral must begin at a time 
T=0 + . The MAM and the FDM produce the same results for a step forcing func- 
tion because for T>0, Q(T) is constant and all its derivatives vanish (see for 
example eq. (16)). As shown in figure 6, the FDM and MAM methods are more accu- 
rate than the MDM and require several fewer modes than the MOM for the same 
degree of accuracy. Figure 7 is a plot of moment error, using the first 25 
modes, as a function of time. Over the time range considered (T=0.001 to 0.01) 
the moment errors of FDM and MAM methods are about one-half the magnitude of the 
error using the MDM method. Also, as time increases the error associated with 
each method for a given number of modes used tends to decrease. 

For a discontinuous forcing function, such as a unit step, y(T) at T=0, the 
modal method exactly predicts a zero response at time T=0. The MAM and FDM 
methods, however, require a summation of all the modes to exactly predict a zero 
response. Therefore, the MDM will produce qualitatively better results for 
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times near T=0 or close to discontinuities. To calculate the transient response 
at very small times a large number of modes is necessary. The MDM using two 
modes, qualitatively approximates the displacement distribution with the same 
accuracy as the MAM and FDM methods using 17 modes for a unit step loading at 
T=.0002 and c-j = .05 (see fig. 8). Figure 9 shows that the error norm associ- 
ated with the MAM and FDM methods is exceedingly large when fewer than about 
five modes are used to approximate the response. If a sufficient number of 
modes are used to predict the displacement distribution accurately (m > 20) the 
FDM and MAM methods appear to give better results. 

For cases where the forcing function or its derivatives are discontinuous at 
some point other than T=0, the appropriate jump conditions must be included to 
use the FDM correctly. The jump conditions are necessary because the integrat- 
ion by parts requires the functions and their derivatives to be continuous. 

This complicates the implementation of the FDM (integrating-by-parts four times) 
for cases where the forcing function is not continuous at least up to its third 
derivative. The FDM can be applied up to the point of discontinuity and then 
re-started using the proper initial conditions. 

Simply-Supported Multispan Beam (10 Spans) 

The second problem studied is a uniform multispan beam (10 equal length 
spans) subject to two loadi ng distr ibutions and one forcing function. The beam 

has a nominal frequency wqVeI/p AL 4 and time, t, is normalized such that 
T = u> 0 t. The first 30 natural frequencies of the beam are normalized and 
listed in table 2. An analytical solution for the mode shapes and frequencies 
of multi span beams was obtained using equations from reference 8. This problem 
was selected because the frequencies are closely spaced (in groups equal to the 
number of spans (10 in this example)) and the chances of a neglected higher mode 
having a considerable effect on the response is increased. 

Case 1: Uniformly distributed load with Q(T)=1000(T 4 -T 5 ) ♦- Figure 10 is a plot 

of the moment distribution, normalized to the maximum value of M, of the multi - 
span beam at T=1.2 and c-j=0.05. The FDM converges using only one mode whereas 
the MDM and MAM methods require 30 and 10 modes, respectively for e<.01 . The 
moment error norm as a function of the number of modes for T=1.2 is shown in 
figure 11. The first nine modes are nearly orthogonal to the uniform load 

distribution, hence, the modal load { 4 >}^{Q(t)} is negligible and has a 

negligible effect on the response (fig. 11). The 10 th and 30 th modes, how- 
ever, have an effect on the solution as shown in the figure. The effect of 
these higher modes, however, is taken into account to some degree, by the 
pseudo- static response (note the MAM curve for m<10) and to a greater degree by 
the higher-order approximation of these neglected modes used in the FDM. Hence, 
the FDM using only one mode gives a more accurate response than the MDM using 49 
modes or the MAM using 9 modes. The moment error norm as a function of time for 
each method (using 10 modes) is shown in figure 12. The accuracy of the FDM is 
at least two orders of magnitude greater than the MDM and one order of magnitude 
greater than the MAM. The MAM and MDM methods are identical at T=1.0 as expect- 
ed because Q=0. At T=1.0, the error associated with the MDM decreases an order 
of magnitude and that associated with the MAM increases an order of magnitude as 
shown. 

Case 2: Two concentrated loads equally spaced about the center of the first 

span with Q(T)=1000(T H -T :> ) .- For this loading distribution, the solution does 
not converge in a step like manner as in figure 11 but does so gradually as 
shown in figure 13. This occurs because the loading distribution is not nearly 


orthogonal to many mode shapes and hence the associated modal load ({♦}{{Q(T)}) 

Is not negligible as In the previous case and these modes contribute to the 
response. Once again, the FDM converges more rapidly than either of the other 
methods and Is at least an order of magnitude more accurate. The accuracy of 
the FDM tends to Increase as T increases (fig. 14). 


Comparison of Convergence for a Sinusoidal Forcing Function 


In this section It Is shown that, for the case of a sinusoidal forcing 
function, the higher-order modal methods such as FDM converge faster than the 
lower-order methods provided that the modal summation encompasses the forcing 
function frequency. For a sinusoidal forcing function ({Q 0 }sinyt) with zero 
initial conditions and no damping, the response using the mode-displacement 
method (MDM) is: 


{U(t)} * { 


2 ){*}{{Q 0 }sinYt 

10 1 

(1 - r^) _ 


{♦} i (*-3) {♦}{ {Q^sinu^t 
"1 


where rf*Y/u>i. The mode-acceleration method (MAM) form of the response can 
be written as: 


{U(t)l . [K] _1 {Q 0 ls1nrt + " ( 

1=1 



Comparing the two equations, the last term in the modal summation is exactly the 
same, the first term in the modal summation is different and will converge 
faster (proportional to r-j 2 ) using the MAM provided that rj<l oraj>Y. 

This is not unreasonable since the frequencies used in the modal summation must 
usually encompass the forcing frequency to provide an accurate transient 
response. The FDM method, using four integration by parts (eq. 16), results in 
the following expression: 


{U(t)} = [Kf^qjsfnrt + [l]' 1 WW' 1 t ! {ll .(slnyt * " ( , 

i*l (1 - r i ) 


{♦} 1 {♦}{ {Q 0 J s<rto i t 

_ H ) 

(1 - r, 2 ) 

The first term in the modal expansion of the FDM converges proportional to 
n 4 or (y/ii)f ) 4 and so for rj<l this term converges much faster than the 
MDM or MAM methods. The last term, however, in the modal expansion remains 
the same for all of the methods. 
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CONCLUSIONS 

The present study evaluates an improved modal method proposed by Leung for 
transient structural analysis and generalizes it to damped systems. The improv- 
ed method, entitled the force- derivative method, is a higher-order method rela- 
tive to the mode-displacement or mode- acceleration methods and is based on 
repeated integrations by parts of the convolution- integral form of the response. 
The repeated integration- by- parts produces terms which give successively better 
approximations to the contribution of the higher modes which are neglected in 
the modal summation. It is shown that the procedure used for deriving force- 
derivative method can be used to derive even higher-order methods and to explain 
the lower-order mehtods such as the mode- displacement and mode-acceleration 
methods. Comparisons are made of the force-derivative, mode-acceleration and 
mode-displacement methods for several numerical example problems for various 
times, levels of damping, and forcing functions. The forcing functions were 
chosen to study the effects of both continuous and discontinuous forcing func- 
tions (with respect to time) and to study what happens when the force or one of 
its derivatives vanishes at some point in time. The example problems include a 
tip-loaded cantilevered beam and a simply-supported multispan beam (10 equal- 
length spans). 

In general, the force- derivative method was found to be more accurate than 
either the mode-displacement or mode-acceleration methods and results in a con- 
verged solution in fewer modes. In addition, for problems in which there are a 
large number of closely-spaced frequencies (e.g., large truss- type space struc- 
tures and multi span beams) the force derivative method is very effective in rep- 
resenting the effect of the important, but otherwise neglected, higher modes. 
Results of a sinusoidal forcing function Indicate that the higher-order modal 
methods will generally converge faster (proportional to the ratio of the forcing 
frequency to the natural frequency raised to some power) providing the natural 
frequencies are higher than the forcing frequency. This is to be expected since 
accuracy suggests that the modal summation Include frequencies which encompass 
the dominant forcing function frequency components. Convergence of the moments 
and shears was found to be slower than convergence of displacements, as expect- 
ed, since the former quantities are functions of the spatial derivatives of the 
displacements. At response times close to discontinuities in the forcing func- 
tion and/or its derivatives, the mode-displacement method qualitatively gives 
better results using fewer modes than the mode- acceleration or force-derivative 
methods. Implementation of the force-derivative method at the discontinuity 
requires the inclusion of the necessary jump conditions. It was also found that 
Increasing modal damping levels does not always increase the convergence rate of 
the mode-acceleration method. The relative importance of damping on the conver- 
gence rate of the modal solution is shown to be complex and dependent upon sev- 
eral factors including the existence of the derivatives of the forcing function. 
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TABLE 1 - NATURAL FREQUENCIES OF A CANTILEVERED BEAM 


(w 


0 



) 


MODE 

NUMBER 

NORMALIZED NATURAL FREQ. 
“ /w o 

1 

3.51601 

2 

22.0345 

3 

61.6972 

4 

120.902 

5 

199.860 

6 

298.556 

7 

416.991 

8 

555.165 

9 

713.079 

10 

890.732 

11 

1088.12 

12 

1305.26 

13 

1542.13 

14 

1798.74 

15 

2075.08 

16 

2371.17 

17 

2687.00 

18 

3022.57 

19 

3377.87 

20 

3752.92 

21 

4147.70 

22 

4562.22 

23 

4996.49 

24 

5450.49 

25 

5924.23 

26 

6417.71 

27 

6930.93 

28 

7463.89 

29 

8016.59 

30 

8589.02 




TABLE 2 


NATURAL FREQUENCIES OF A SIMPLY-SUPPORTED 
MULTISPAN BEAM (10 SPANS) 



MODE 

NUMBER 

NORMALIZED NATURAL FREQ. 
“ /(0 o 

1 

9.86960 

2 

10.1479 

3 

10.9435 

4 

12.1691 

5 

13.6959 

6 

15.4212 

7 

17.2490 

8 

19.0691 

9 

20.7508 

10 

21.9115 

11 

39.4784 

12 

40.0728 

13 

41.7202 

14 

44.1072 

15 

46.9043 

16 

49.9649 

17 

53.1222 

18 

56.2347 

19 

58.9525 

20 

60.9470 

21 

88.8264 

22 

89.7765 

23 

92.2339 

24 

95.7363 

25 

99.8682 

26 

104.248 

27 

108.879 

28 

113.289 

29 

117.133 

30 

119.972 
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a) Displacement distribution 

Figure 1.- Representation of displacement and moment errors In a cantilevered 
beam with tip loading using error norm e (Q(T) * 1000(T 4 -T 5 ), 

T - 0.4 and * 0.05). 
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Figure 2.- Comparison of moment errors of a cantilevered beam subject to a tip 
load Q(T) * 1000(T I, -T 5 ) (T=1.2 and c-j = 0.05) using three 
different modal methods. 
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Figure 4.- Comparison of displacement, moment, and shear errors of a 
cantilevered beam using three different modal methods (ci 
and tip load Q(T) = 1000(T**-T 5 )). 
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Figure 6.- Comparison of moment errors of cantilevered beam subject to a unit 
step tip load Q(T) = y(T) (T = 0.01) and = 0.05) using the 

different modal methods. 
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Figure 7.- Variation of moment errors as a function of time , assuming 25 modes 
are used in the modal summation, for a cantilevered beam with a unit 
step tip load Q(T) = W (T) ( Ci = 0.05). 
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Figure 8.- Displacement distribution In a uniform cantilevered beam with a unit 
step tip loading at time T = 0.0002 (Q(T) = p(T) and 
= 0.051. 
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Figure 10.- Comparison of normalized moment distributions along a 

simply-supported multispan beam (ten equally-spaced spans) and a 
uniformly distributed load Q(T) = 1000 (T 4 -T 5 ) (T = 1.2 and c-j = 
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Figure 11.- Comparison of moment errors of a simply -supported multispan beam 
(ten equally -spaced spans) subject to a uniformly distributed load 
Q(T) = 1000 (r*-T 5 ) (T « 1.2 and cj = 0.05). 
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Figure 12.- Variations of moment errors as a function of time, using ten modes 
in the modal summation, for a simply -supported multi span beam (ten 
equally-spaced spans) subject to a uniformly-distributed load 
Q(T) = 1000(T 4 -TS) (ti * 0.05). 
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Figure 13.- Comparison of moment errors of a simply-supported multispan beam 

(ten equally-spaced spans) subject to two concentrated loads spaced 
about the center of the first span and varying as 
Q(T) = 1000(T 4 -T 5 ) (T = 0.4 and Ci = 0.05). 
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Figure 14.- Variation of moment errors as a function of time, using 10 modes in 
the modal summation, for simply-supported multispan beam (ten 
equally-spaced spans) subject to two concentrated loads spaced 
about the center of the first span and varying as 
Q(T) * lOOOCr 4 -! 5 ) ( C1 = 0.05). 
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